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ABSTRACT 

We present constraints on cosmological and astrophysical parameters from high-resolution mi- 
crowave background maps at 148 GHz and 218 GHz made by the Atacama Cosmology Telescope 
(ACT) in three seasons of observations from 2008 to 2010. A model of primary cosmological and 
secondary foreground parameters is fit to the map power spectra and lensing deflection power spec- 
trum, including contributions from both the thermal Sunyacv-Zcldovich (tSZ) effect and the kinematic 
Sunyacv-Zoldovich (kSZ) effect, Poisson and correlated anisotropy from unresolved infrared sources, 
radio sources, and the correlation between the tSZ effect and infrared sources. The power (^Ci/in of 
the thermal SZ power spectrum at 148 GHz is measured to be 3.4 ± 1.4 /iK^ ad, £ = 3000, while the 
corresponding amplitude of the kinematic SZ power spectrum has a 95% confidence level upper limit 
of 8.6 ^K"^. Combining ACT power spectra with the WMAP 7- year temperature and polarization 
power spectra, we find excellent consistency with the LCDM model. Wc constrain the number of 
cifFective relativistic degrees of freedom in the early universe to be iVeff = 2.79 ± 0.56, in agreement 
with the canonical value of A^eff = 3.046 for three massless neutrinos. We constrain the sum of the 
neutrino masses to be Sm^ < 0.39 eV at 95% confidence when combining ACT and WMAP 7-year 
data with BAO and Hubble constant measurements. We constrain the amount of primordial helium to 
be 1^ = 0.225 ± 0.034, and measure no variation in the fine structure constant a since recombination, 
with a/ao = 1.004 ± 0.005. We also find no evidence for any running of the scalar spectral index, 
dus/dlnk = -0.004 ±0.012. 

Subject headings: Microwave Telescopes, CMB Observations 
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1. INTRODUCTION 



Studies of the cosmic microwave backgromid (CMB) 
have dramatically progressed over the past two decades 
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2003; Johnson ct al. 2007; Chiang ct al. 2010). The cur- 
rent ACDM cosmological model provides an excellent fit 
to the CMB data across a wide range of angular scales, 
and is supported by c omplementary observations of larg e 
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Collaboration et al.|2012 1 . As this paper was being final- 
ized the SFl' collaboration releas ed a new set of papers 
( [Story et al.l2012l|Hou et al. | 201 2 ]) and the WMAP team 
released its final 9-year results ( [Bennett et al.|2012 Hin- 
shaw et al. 2012). Our analysis does not incorporate 
these results although in a few places we make direct 
comparisons. We anticipate the release of the cosmolog- 
ical results from Planck soon. 

The Atacama Cosmology Telescope (ACT) comple- 
ments measurements from WMAP by observing from £ ~ 
300 to £ = 10000. This widens the range of data avail- 
able to constrain both cosmological paramete rs through 
the Silk damping tail of the primary CMB ( [Silk[[1968 1 
and the residual power from secondary sources between 
us and the surface of last scattering. These sources in- 
clude galaxy clusters, which are detectable at microwave 
iiiuK-eii ei ai ^uu» ivessiei eu ai ^uua frequencies through the Sunyacv-Zcl'dovich (SZ) effect 
20l !i',T a"axrcj l^ £ersm^ gel'dovich fc Sunyaev 1969, Suny aev fc ZePdovich 1970 ) . 

Hionnn, A/ r ^i l omnllp^r,^ ^i l oni n l Ihe thermal SZ (tSZj ettect describes the spectral distor- 
tion due to the inverse-Compton scattering of CMB pho- 
tons to higher frequencies by the hot gas in clusters, while 
the kinematic SZ (kSZ) effect measures the correspond- 
ing temperature shift due to the bulk peculiar motion of 
the clusters. While these effects produce a diffuse signal 
at small scales from unresolved clusters, their influence 
has recently been detected directly through the cross- 
correlati on of the ACT t e mper ature maps with other 
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Tinker et al.|201^ ) and observations of gravitation al lens 
ing (e.g.,[Massey et al.|l2007||Fu et al..2008, .Schrabback 
eraL][20lOI [Suyu et al. 2010; Hcyn Tans et al.|[i^01ij| [Kil- 
binger et al.[[2012 ). While the CMB on angular scales of 



greater than 0.3" has been definitively measured b y the 
Wilkinson Micr owave Anisotropy Probe (WMAP, Ben- 
nett et al.|2012 ), a wealth of information in the CMB on 
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tracers (Hand et al. 2011 [2012|. In addition to emis 
sion from clusters via the SZ effect, radio galaxies and 
dusty star-forming galaxies also contribute to the power 
on small scales, and indeed dominate the cosmological 
signal for multipoles £ > 3000 and frequencies i/ > 90 
GHz. Gravitational lensing by structures along the line 
of sight also generates a microwave background signal, 
and distorts the primordial CMB. 

This paper forms part of a set presenting the 3-year 
analysis of the ACT data; the ACT te mperature and de - 
flecti on power spectr a are p resented in Das et al. (2013), 
while Dunkley et al. (|2013[) presen ts the likelihood used 
in this analysis. [Hass elheld et al.| ((2013|) presents a cat- 
alogue of SZ-detected clustcrs^ro^l tire ACT data, and 
interprets them. This paper contains the parameter es- 
timation of both primary (cosmological) and secondary 
(foreground) parameters. We outline the data used in 
this analysis in Section [2| and describe the methodol- 
ogy and likelihood in Section [3] We present constraints 
on primary cosmological parameters in Section l4| and on 
secondary parameters in Section [5| We conclude in Sec- 
tion [6] after which we provide an appendix of analysis 
tests. 

2. DATA 

This paper presents results from a combination of ob- 
servations at two frequencies, 148 GHz and 218 GHz, of 
multiple fields taken over three years. The southern fields 
(ACT-S) varied over different seasons, with the 2008 sea- 
son containing a 292 deg^ patch and the 2009 and 2010 
seasons focused on a smaller 146 deg^ footprint. Equa- 
torial data (ACT-E) were taken only over the 2009 and 
2010 seasons; we use a 300 deg^ patch for the 3-year 
analysis. We follow a similar procedure for going from 
maps to the temp erature power spectra as was used in 
Das et al. (2011b), usi ng the power spec trum estimation 
procedure presented in Das et al. ( 2009 1 . 
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Fig . 1. — Data used in this cosmological analysis. The data from the WMAP 7-year data release ( [Larson et al.||201l ! 'Komats u et al.| 
|2011[ | are combined with the ACT data. In this figure we show a weighted co-added spectrum from the equatorial and southern patcUes at 
148 GHz. The full ACT likelihood, however, considers independent southern and equatorial spectra for both 148 GHz and 218 GHz (given 
in Table 3 of Das et al. 2013 , and online on LAMBDA) and their cross-frequency spectrum, which are shown in Figure pi The solid line 
indicates the best-ht cosmological model including foreground emission, while the dashed line shows the best-fi t primordiaFCMB spectrum. 
This binning was selected for the common analysis of the equatorial and southern data. See|Das et al.|(|2013^ for alternative binnings. 



The data and map -making procedure are described in 
Diinner et al. (2013), the powe r spectrun i meth od and 
systematic tests are prese nted in| Das et al. (2013). Using 
the spectrum presented in]Das et al. (2013|, we have con- 
structed two hkcUhoods; these likelihoods are presented 
in Dunkley et al.^ (2013) . The multi- frequency likehhood 
parameterizes the foreground emission using additional 
parameters, which we hst in Section [3j while the CMB- 
only likelihood marginalizes over these foregrounds. For 
the analysis presented here we use the multi-frequency 
likelihood; however we show in Appendix [D] that the two 
likelihoods give equivalent results. 

The data used in this analysis are the multi-frequency 
temperature spectra estimated from the ACT maps with 
combined spectra from both ACT-S and ACT-E in clud- 
ing their c ovariance. The spectra are presented in |Das| 
et al. (|2013). For the results presented in Figure ji] we use 
the rnarginalized CMB-only likelihood. In addition, we 
include the measurement of the p ower spectrum o f the 
lensing deflection angle, P-j^Ci''^ ( jPas et al.||2013j) 
The lensing spectrum is estimated from the AUl' tem- 



perature maps usi ng an optimal quadratic estimator ( Hu 
fc Okamoto 2002). Only data from the ACT equatorial 
patches are used to measure the deflection power. The 
covariance between the lensing power spectrum and the 
temperature power spectrum is small. When adding in 



the deflection data, we will use the abbreviation 'ACT- 
Defl'. 

The maps are cross-correlated with WMAP7 maps in 
order to obtain a calibration factor in multipole space. 
The cross-correlation calibration method is described in 



Hajian et al. ( 2012 ); det ails of the calibrat ion of the ACT 
3-year data are given in Das et al. (2013). 



2.1. Beam and calibration errors 

Understanding the beam profiles is essential for inter- 
preting the high-.^ aspects of the power spectrum. At 
e = 3000 and 148 GHz, the window function is w 0.5 
its value at £ = 200. The beams ar e estimated indepen- 
dently for each array and season (Diinner et al. 20131 
from observations of Uranus. The beams vary slightly 
with season due to changes in the telescope focus. We 
include a contribution in the likelihood to the full covari- 
ance matrix from the covariance of the beams for each 
season. The beam error includes contributions from the 
uncertainty in the pointing variation of the telescope. 
At £ = 700, the pivot point for the beam and calibra- 
tion uncertainties, the effective calibration error is 2% 
for the 148 GHz maps and 2.6% for the 218 GHz maps. 
Different seasons have somewhat different calibration un- 
certainties and so these numbers should be considered as 
representative of the effective combined calibration. The 
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calibration uncertainties of ACT-E and ACT-S are sta- 
tistically independent. 

In the cosmological analysis, we apply a calibration 
prior for the 148 GHz spectra obtained from the 
WMAP-ACT cross-correlation calibration procedure 
described above. In addition, we allow for a small 
error in the overall calibration of the spectrum by 
marginalizing over calibration factors for the south and 
equatorial spectra, at both 148 GHz and 218 GHz. 
The extra calibration allows for the overall spectra to 
adjust themselves at the 1% level and has a negligible 
effect on the cosmological parameters, as is discussed 
in Appendix [C] Similarly, in the same appendix, we 
test for the dependence of the cosmological parameters 
on the beam assumed uncertainty, and find that beam 
error has a negligible effect on the parameters of interest. 



2.2. Additional data 
We use temperature and polarization data fr om the 



seven year data release of the W MAP satellite ( Larson 
et al.|20lT Komatsu et al.|2011 1 in addition to the mea- 



surements oFTEie'lmcrowiiTC from ACT. We 

include measurements of the Baryon Acoustic Oscilla- 
tions (BAO) from the Six-degree Field Galaxy Redshift 
Survey (6dFGRS, 'Beutler et al.||2011[) and the Sloan 
Digital Sky Surv ey Data Releases 7~[ BDSS DR7, [Per- 
cival et al.||2010[ ) and 9 (SDSS DR9, [Anderson etaT 
2U12|), measured at redshifts: z = 0.10t),U.35 and U.57. 



in addition, we supplement our data with the latest 



measurement of the Hubble co: 
2.4 km s"^ Mpc"^ (Riess et al. 



man et al. 



Slant of Ho = 78.8 + 



(2012) and Ho 



2011j), although [Freed- 
74.3 ±2.1 km s-i Mpc 



In 



some cases we also include a prior on ag f rom skewness 



meas urements of the tSZ effect from ACT (Wilson et al. 
20T1. 



Unless explicitly specified, the ACT 3-year data are 
combined with WMAP7 data. For some model con- 
straints, such as those on N^g and the secondary param- 
eters, we also include the 'low-£' and 'h igh-£' spectrum 
measurements from SPT ( Keisler et al.|p0 11; Rcichardt 
et al.| 2012), following the prescription in pJunkley et aTj 
(fmgp We show the ACT and WMAP7 data in Figure[lj 
'i'he best-fit model for the various frequency components 
is shown in Figure [2j 

3. METHODOLOGY 

We use Markov Chain Monte Carlo (MCMC) meth- 
ods to determine parameters associated with a variety 
of models. The basic cosmological ACDM model con- 
sists of 6 parameters describing a flat universe. These 
include the physical baryon density fl^h'^ (where h = 
iJo/100 km s~^ Mpc~^ is the dimensionless Hubble pa- 
rameter), cold dark matter (CDM) density ftcfi^, and 
Oa = ra/dA, the ratio of the acoustic horizon to the an- 
gular diameter distance at decoupling. This parameter 
is sensitive to the dark energy den sity, S7a, but less de - 
generate with other parameters ( jKosowsky et al. 2002). 
The value of fl^ is then a derived parameter. We assume 
the primordial perturbations to be scalar, adiabatic, and 
Gaussian and parametrize them via a spectral tilt rig, 
and amplitude A^, defined at pivot scale ko = 0.002 
Mpc~^. We assume that the universe transitioned from 



a neutral to an ionized state over a small redshift range, 
Az = 0.5, with optical depth r. The reionization history 
of the universe can be probed by small-scale CMB mea- 
sureme nts through the impact of reionization on the kSZ 
effect (Ostriker 
[Knox ct al 



Vishniac 1986 



Gruzinov fc Hu||1998 
, although care must 



1998; /ahn et al.,,20I 

be taken to allow tor correlations between the tSZ ef- 
fect and the microwave emission from unresolved du sty 

[Addison et al.||2012b[ ). 



galaxies (Mesinger et al. 



We express the ACDiV 



2012 



set o: 



parameters as 



{^76/l^^]c/i^6'A,A^,n,,T}. 



(1) 



In this ACDM model, the number of effective relativistic 
degrees of freedom is assumed to be N^tf = 3.046, with 
the abundance of primordial helium fixed at Yp = 0.24. 
The likelihood used i n the ACT analysis is described in 
(2013), which we briefly summarize here, 
of secondary emission to the ACT multi- 



Dunkley et al. 
We ht a mode. 



frequency power spectra that includes an additional nine 
parameters when considering ACT data in combination 
with WMAP7 data. For the 148 GHz data we use modes 
500 < £ < 10000, while for the 218 GHz data we restrict 
ourselves to 1500 < £ < 10000. The theoretical spectrum 
for frequency bands i and j is 



th.ij 



CMB 



(2) 



where Bg = £{£ + l)Q/27r and C^^^ is the lensed pri- 
mary CMB power spectrum. The secondary spectra 
components are modeled as 



gscc,ij^gtSZ,ij_^gkSZ,ij 



-B 



CIB,ij 



B 



tSZ-CIB,ij 



B- 



rad,ij 



B 



Gahij 



(3) 



with contributions from the tSZ and kSZ effects, GIB 
sources, the cross-correlation between the tSZ and GIB 
signals (tSZ-CIB), radio galaxies (rad), and residual 
Galactic dust (Gal). In addition to the six primary pa- 
rameters, we add the following nine parameters 



{fltSZ, flkSZ, flp, flc, 0,gs,age, Pc, Os, C}- 



(4) 



The parameter otsz parameterizes the amplitude of the 
tSZ power; Oksz the kSZ amplitude; Op and Oc the Pois- 
son and clustered Cosmic Infrared Background (GIB) 
power and agsjflge model the residual Galactic dust 
anisotropy in the southern and equatorial survey regions. 
All the parameter amplitudes are dimensionless, and are 
defined for a template spectrum normalized to 1 /iK^ at 
£o = 3000, and frequency v = 150 GHz. The frequency 
dependence of the correlated and Poisson GIB power is 
given in flux density units by the product of modified 
blackbodies with effectiv e temperature 9.7 K a nd emis- 
sivity index /3c, following Addison et al. (|2012c )_^and de- 
scribed in equations (8) and (9) of Dunkley et al. (2013). 
The radio source power has an amplitude and a spec- 
tral index as fixed to —0.5, based on the assumption 
that the index obtained fo r brighter sources from ACT 



and SPT so urce catalogs (Vieira et al. 2010| [Marriage 



et al.||2011a ) holds for fainter sources. 

'I'he clustered and Poisson (both GIB and radio) tem- 
plates vary with scale as £{£ + l)Cf/27r cx £^-^ and cx £^, 
respectively. We allow for a correlation between the tSZ 
effect and GIB sources, wit h scale dependence giv en by 
the template calculated by Addison et al. (2012b), and 
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a frequency-independent corr elation coefBcient, w hich 



is defined in Equation (11) of Dunkley et aL| (|2013[), and 
is restricted to he in the range U < ^ < (J.2. I'ne parame- 
ters in our foreground model are summarized in Table [l] 

The spatially variable Galactic emission has been 
masked, leaving only a small residual component. To de- 
termine the amplitude and spectrum of the resid ual com- 
ponent we cross correlate with the IRIS map ( Miville- 
Deschenes fc Lagache|2005 1 as described in detail in Das 



et al.| puia p. We then marginalize over this dust com 
ponent m the likelihood with separate amplitudes in the 
equatorial and southern regions. There is roughly double 
the amount of dust in the equatorial region than in the 
south. In addition, there are some bright clouds. Investi- 
gations were performed to identify possible residual dusty 
clouds that were missed by our treatment, but no clear 
evidence for them was found. The likelihood marginal- 
ization described in Dunkley et al. (2013) was found to 
be the most general and parsimonious treatment. 

We mask all sources above a detection threshold of 15 
m Jy. The source d etection algorithm is described in |Mar-| 
riageetal. (201 lb I. The 148 and 218 GHz source s amples 
in the somhern niap will be presented in Marsdc n et al.| 
(2013), and source catalogs for the full data set will be 



presented in Gralla et al. (2013). Sources are masked to 



a flux level ot 6 . 4 mJy in the SPT analysis presented in 



(2011) 'low-^' spectrum. 



9.2 fxK^ from the Keisler et al. 

For the parameter analysis, we use the publicly avail- 
able CosmoMC code, which includes version 1.5 of the 



Recfast code (fScagcr et al. 1999, 2000; Wong et al. 


2008 


Switzer & Hif ata, ,2008 


,Ali-Haimoud & Hirata |2011 


Chluba & Thomas (2011 


)• 



lensed cross spectra: 148x148 GHz, 148x218 GHz and 
218x218 GHz made from the ACT-S and ACT-E fields 
described in Section 2] In addition, we use the measure- 
ment of the lens ing deflection field by ACT, presented in 



( Das et al.|2013 ). We investigate four types of fits in this 



analysis: 



• We apply the full ACT likelihood to the data in 
combination with WMAP7 data and obtain con- 
straints on both the primary and secondary param- 
eters, with a full 15, 16 or 17 parameter model (and 
four additional calibration nuisance parameters). 

• We estimate Ci bandpowers marginalized over the 
secondary foregrounds, and obtain constraints on 
primary parameters based on these marginalized 
bandpowers. 

• We combine ACT and SPT to check for consistency 
between these small-scale experiments. 



|Reichardt et al. ( 2012 ) , resulting in a source amplitude at 
150 GHz ot a ' = 1.3 ji 0.2. Sources are masked to a level 
of 50 mJy in Keisler et al. (2011). The estimated differ- 
ence in unresolved radio source power is 9.2 /iK^. Hence 
we include a separate amplitude for radio emission mea- 
sured by the South Pole telescope, a'^, when including 
the SPT data in our analysis, whereas the other param- 
eters in the secondary model, apart from the Galactic 
dust parameters, are common to both data sets. Thus 
we first subtract an amplitude of radi o Poisson power of 
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Fig. 2. — The ACT multi-frequency data. Each soUd line indi- 
cates the best-fit total model for each cross spectrum, while the 
dashed lines show the radio power; the dot-dashed line indicates 
the power from CIB sources (both Poisson and correlated CIB); 
the gray lines show the total SZ power. The solid and dotted gray 
lines are for those spectra including 218 GHz and consist of only 
the kSZ contribution, while the dashed gray line is for the 148x148 
GHz spectrum. The dotted lines show the (negative of) the (SZ- 
CIB) correlation. In all cases the models are for the best fit in a 
ACDM scenario. 

• We consider the ACT data alone without WMAP7 
data, while fixing (or placing priors on) the spectral 
index ris and the optical depth r. 

The models obtained from fitting only the ACT-S and 
ACT-E spectra are consistent with the models fit to 
the combined data, the best-fit spectra from each region 
agreeing to within 4%. Moreover, the best-fit theoreti- 
cal spectrum from the current 3-year ACT data agrees 
with the spectrum derived from the 1-year ACT data to 
the 1% level. We discuss the consistency of the ACT 
spectrum in Appendix \K\ 

4. CONSTRAINTS ON PRIMORDIAL PARAMETERS 

The ACDM model continues to fit the ACT data 
well, when combined with the independent WMAP7 
data. Figure Is] illustrates the constraints on the ACDM 
model with the additional secondary parameters for the 
ACT-I-WMAP7 data combination. In addition, we plot 
the constraints from the WMAP7 power spectrum alone. 
ACT extends the angular range measured by WMAP, 
but the parameters from the joint fit are consistent with 
those from WMAP alone. In addition, the plot shows 
that the six parameters are robust to the presence of low 
levels of foreground emission that can be identified and 
extracted by ACT because of its higher resolution and 
different frequency coverage. 

We start by constraining the parameters in the ACDM 
model. Our constraint on the scalar spectral index is 
= 0.971 ±0.009, using ACT data in combination with 
WMAP7 data, BAG and Hq mea surem ents. We discuss 
the constraints on rig in Section |4.5| In addition, we 
improve the constraints on the baryon density fli,h^ using 
only CMB data, to ^bh^ = 0.0225±0.00047, which is due 
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TABLE 1 

Secondary parameters in the ACT 3- year foreground model. 



Symbol 


Description 




Thermal Sunyaev-Zel'dovicli (SZ) power amplitude. 




Kinematic SZ power amplitude. 




Poisson Cosmic Infrared Background (CIB) power amplitude. 


ac^ 


Clustered CIB power amplitude. 




Residual galactic emission amplitude for the ACT-S spectrum. 




Residual galactic emission amplitude for the ACT-E spectrum. 




Emissivity index of the clustered CIB power. 




Radio Poisson power amplitude. 


e 


tSZ-CIB correlation amplitude. 



^Amplitudes relative to templates normalized to £{£ + l)Ci/2n = IfiK^ at £q = 3000 and 150 GHz. See [Punkley et al.] j2013[ l for a 
complete description. 




0.03 0.09 0.15 456789 10 2 3 4 13579 




WMAP7 (Recfa.stvl.5) 

ACT+WMAP7 

ACT+WMAP7+BAO+H0 



1.7 2.0 2.3 2.6 0.0 0.2 0.4 0.6 0.8 1.0 0.2 0.6 1.0 1.4 



Fig. 3. — Parameter constraints on the ACDM model for the combined ACT and WMAP 7-year power spectra. For each parameter, the 
full likelihood is marginalized over the other parameters. Dotted curves for the primary parameters are for WMAP 7 data only. The solid 
curves are fo r ACT plus W MAP7. and the dashed lines add priors from measurements of the baryon acoustic peak ( [Anderson et al.|2012| 
|Perc_ival et a l. 2010 Bcutler et al.|2011[ | and a direct measurement of Hq (Riess et al. 2011). The constraints are sunimanzed m Tables [^p| 
and[b] bmce the § parameter is unconstrained and set only by the prior, < ^ < U.2, it is not shown in the tables. 



to the fact that the ACT spectrum places tight hmits on 
the positions of the higher order peaks below £ « 3000. 
The models are summarized in Tables [2] and |4] for the 
ACT data in combination with WMAP7. Table [5] shows 
the constraints when adding BAO and Hq measurements 
to the ACT power spectrum data. In Appendix ID] we 
show that the ACDM parameters derived from the CMB- 
only likelihood agree with the full likelihood to within 
0.1a. 



The greatest power of ACT comes when quantify- 
ing models beyond the standard cosmological model be- 
cause the temperature power spectrum contains little ad- 
ditional information on the simple ACDM parameters 
at angular sca les smaller than the th ird acoustic peak, 
i « 800 (e.g., [Kosowsky et al.||2002[ ). The damping of 
the hi gher-order acoustic p eaks relative to the baseline 



model ( Komatsu et al. 2009 ) provides constraints on a va- 
riety of non-standard models. The ACT 1-year spectrum 
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Fig. 4. — Constraints on the ACDM cosmological model using the likehhood of marginahzed CMB bandpowers using ACT+ACTDefl data, 
compared to the constraints from ACT+ACTDefl+WMAP7, where ACTDefl is the reconstructed deflection power spectrum from ACT 
data, and WMAP7 alone. The top six panels are fltted parameters, while the bottom two rows are derived parameters. For the ACT-only 
parameters, two cases are shown: one where we fix the optical depth r, and another where the scalar spectral index Us is fixed. Fixing ris 
tightens the constraint on the amplitude of scalar perturbations (the ACT+ACTDefi, ris fixed curve lies on the ACT+ACTDefi+WMAP7 
curve), but does not cause significant shifts in the other primary parameters, relative to leaving ris free. The independent CMB data sets, 
ACT and WMAP7, are consistent with the ACDM model. 



data ( Das et al. 2011b ) showed a slight excess of damping 
at small scales relative to the baseline model. Evidence 
for this slight excess is not seen in the ACT 3-year data 
set. More data resulted in a spectrum with smaller error 
bars, providing tighter constraints on parameters such as 
the baryon density, while the best-fit theoretical spectra 
are consistent between the two results at the 1% level. 
The consistency b etween the results pre sented here and 
those presented in Dunkley et al. (2011) is discussed in 
Section |D] of the appendix. 
Various authors have explored possible models which 



Galh et al.||2011| |Calabrese et al.||2011a 


b; Hasenkamp 


2012 


HamannI 2012 Menestrina & Schcrrer 2012 |Foot 


2011 


lAbazajian et al.| |2011| |Menegoni et al.| |2012| 


Farhang et al. 2012p. In this analysis we broaden our 



the damping tail of the ACT data in this context. 



4.1. ACT data alone 

In Figure |4j we compare the constraints on the 
ACDM mode" in the case of ACT data alone com- 
pared to the WMAP7 constraints, and the combination 
of ACT-I-WMAP7. This provides an important cross- 
check of the consistency of the two data sets. A key 
parameter which is primarily constrained with polariza- 
tion data on the largest scales, such as those probed by 
WMAP, is the optical depth r. Hence, when considering 
the constraints from the ACT data alone, without in- 
cluding WMAP7, we impose a prior on the optical depth, 
T = 0.089 ±0.015. In Figure |4| we show two cases, one 
in which the scalar spectral index is fixed at Us = 0.967 
and one in which rig is allowed to vary freely. Fixing 
rig tightens the bound on the amplitude of fluctuations, 
while other parameters are largely insensitive to the ef- 
fect. The agreement shows that the same model that 
describes the WMAP7 data for i < 1000 independently 
fits the damping tail measurement from ACT o f £ > 5 00. 
The same behavior is observed for Story et al. (2012). 
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Fig. 5. — Marginalized one-dimensional likelihood of the N^ff for ACT data in combination with other probes. Left panel: The improve- 
ment in constraints of the ACT data in combination with WMAP7, relative to WMAP7 data on its own. i?ight panel: The N^g likelihood 
for WMAP7+ACT data and a variety of other probes. The constraints are consistent with A^gfj = 3.046, the value in the model with 
three neutrino species. The slight shifts in the central value can be understood in the context of the changes in the value of the Hubble 
parameter, as illustrated in Figure [6] 



4.2. Effective number of relativistic species 

The standard cosmological model has three neutrino 
species, all of which have negligible mass and contribute 
to Ncs, the effective number of relativistic species at 
recombinationlf3 Precision electroweak measurements 
place tight constraints on the number of light neutrino 
species with standard-rn odel couplings (The ALEPH 
Collaboration et al.|2006 ) through Z production in e+e" 
collisions: 

A^eff = 2.984 ± 0.008. (5) 

Relativistic species (be they neutrinos or another rel- 
ativistic species at early times) change the expansion 
rate of the universe through their energy density and 



the damping tail of the primary CMB spectrum 1 


Bowen 


et al. 


2001 


Bashinsky & Seljak 


2004 


Hou et al. 


2011|). 



than that of photons by a factor 

p,/p^ = (7/8)(4/ll)4/3Ar^ff, 



(6) 



where N^s is 3.046 in the standard ACDM model. 

Extra relativisti c energy density d amps the small-scale 
CMB power - see | Hou et al.| ( [20lTl) f or a concise recent 



review, and the discussions in Hufc Dodolson (2002) 



Hu et all (1 20011); 'Bashinsky fc Seljak (2004); Tcgmark 
( |2005[ );|Lesgourgucs & Pastor (2006); Hanncstad (2010). 

l'igure[5| and Table [2] illustrate the constraints on the 
A'off from~ACT in combination with various probes. Pre- 



nalyses ( iDunkley et al.||20lH [Keisler et aL||2011 1 
,ed a slight excess in the iVpff . i'his preference tor 



vious a: 

suggested a slight excess in the I\es- i'his prete: 
more damping from extra relativistic degrees of freedom 
is no longer present when analyzing the ACT 3-year data, 
in combination with WMAP7 data. The change is con- 

The effective number of relativistic species due to three neu- 
trinos is slightly larger than three even in the standard scenario 
due to heating caused b y the injection of the entropv from the 
e+/e~ annihilation (e.g., Dicus et al. 1982, Rana & Seifert 'l99l' 
Dodels on & Turner 1992 Uolgov & I'ukugita 1992 Hanncstad & 
Madsen 1995j |Uolgov et al.||1997| |Gnedin fc Gnedin||199S( |Lop^ 

* X y y y ) . 



sistent with the improved statistics of the ACT 3-year 
data. We find 



off 



2.79 ±0.56 (WMAP7-f ACT). 



(7) 



The improvement in this value relative to the WMAP- 
only constraints is shown in the top panel of Figure [6j 

The result in this analysis was obtained by imposing 
the consistency relation between the primordial helium 
fraction at Big Bang Nucleosynthesis (BBN) a nd the 
number of e f fective relativistic degrees of freedom (|Trotta 
eTaT] [2004 1, 'Kneller & Steigman '2004; ' SteigmariP ' 
Simha TTSt cigman 2008, see Section [4. 9[); 



Yp = 0.2485 + 0.0016[(273.9O6/i2 - 6) + 100(5* - 1)]; 
5= Vl + (7/43)(Afeff-3). 

(8) 

Hence, in the present analysis, the helium fraction is a 
determined parameter given N^g and the baryon den- 
sity, rather than remaining fixed at the standard value of 
Fp = 0.24. The value presented in IDunk ley et al.| ( |2011 | 
was higher at N^g = 5.3 ± 1.3, as we did not impose this 
relation; imposing this constraint on the previous ACT-S 
data would yield a modified value of A^eff = 4.3± 1.3, con- 
sistent at la with the value of 3.046 expected in standard 
ACDM. 

The value of Acff obtained when using only the ACT- 
S data is closer to the value reported in |Dunkley et al.| 
( |2011| . The ACT-E data prefer a lower value for A^eff, 
leading to a combined resu l t whic h is w la lower than 
presented in Dunkley et al. (20111. 

Including the recent BAO data does not change the 
constraints on the relativistic species, while adding in the 
SPT data shifts the mean value to slightly higher values 
of A^cff , but stiU consistent with the WMAP7+ACT data. 
The highest value for Ngg comes from the addition of the 
BAO and Hubble constant data, yielding 

N^ff = 3.50 ± 0.42 (WMAP7 -I- ACT + BAO + HST). 

(9) 

This is due to the degeneracy between ilch"^ (and there- 
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Fig. 6. — Marginalized two-dimensional (68,95%) contours in 
the N^g — Hq plane. The two parameters are highly correlated, 
with larger values of Hq leading to higher values of the effective 
number of relativistic degrees of freedom. Here BAO refers to the 
combined 6dF and SDSS DR9 BAO measurements, while HST 

r 



is the Ho = 73.8 ± 2.4 km s Mpc measurement from Rioss 
|et al.| |2011[ l. In contrast, the analysis of Chen & Ratra (2011) 
prefers Hq = 65.5 it 5.5 km s"^ Mpc~^, where the 95% error bar 
includes both systematic and statistical errors. When combining 
both BAO and HST measurements, TV^ff is pushed to higher values 
than in either of the individual cases, as N^g tries to reconcile the 
differences in the distance ratio rs/DY (where rs is the co-moving 
sound horizon at recombination and Dy is an effective geometric 
distance) between the two probes. A similar trend is seen for the 
matter density, f2c^^- 



fore the Hubble constant in a flat universe) and the rel- 
ativistic species, shown in Figure [61 In a flat universe, 
higher Hq leads to lower flch^, which increases the power 
on medium to small scales £ ^ 200 (as the radiation driv- 
ing of the acoustic oscillations is reduced) , leading in turn 
to a larger value of Ncs needed to damp power in the 
tail of the spectrum. The mild tension between the Hq 
inferred through BAO distance measurements and the 
Hubble constant measurements leads to a value of Ncs 
which is higher than the constraint when only the Hub- 
ble constant is added to the CMB data. A similar trend 



was s een in the recently released SPT results ( |Hou et al. 



20121. Smaller valu es of the Hubble constant prefer lower 



Chen & Ratra 


2011 


are summarizec 


in 1 



Calabrese et al. 



In addition, the correlation between the scalar spec- 
tral index and the iVeff is shown in Figure [Tj Decreasing 



Fig. 7. — Marginalized two-dimensional (68, 95%) contours in the 
A'eff — ris plane. The two parameters are strongly correlated, with 
larger values of Ub leading to higher values of the effective number 
of relativistic degrees of freedom. This figure is based on the 16 
parameter fit of the standard cosmology and including variations 
in N^ff. 



power at small scales (through increasing Ncs) is com- 
pensated by increasing the scalar spectral index, which 
increases small-scale power. 

Finally, we also consider an additional constraint from 
the recent ACT measurement of the skewness induced 

The tSZ skew- 



2012). 



by the tSZ effect ( [Wilson et al.[ 
ness signal is more sensitive to ag than any other cosmo- 
logical parameter (scaling approximately as (Jg^), allow- 
ing for a tight constraint with few degenera cies. Using 
theoretical calculations similar to those in Wilson et al.l 



(2012), we find that the most significant degeneracy is 
with ^bh? , for which the tSZ skewness scales approxi- 
mate ly as (^bh'^Y' ^ . Thu s, we include the erg constraint 



from Wilson et al.j ( 2012[ ) in the following form 
cr8(rj6/iV0.0226)°-3 = 0.79 ± 0.03, 



(10) 



where 0.30 = 3.3/11.1, and 11.1 is the fiducial scaling 
of the tSZ skewness value with erg. In addition to its 
correlation with n^, the effective number of relativistic 
degrees of freedom is strongly correlated with trg: as N^fi 
increases, so does CTg, hence this prior lowers the effective 
number of relativistic degrees of freedom to A^eff = 2.55± 
0.40. 

The fact that ACT resolves the higher order peaks of 
the CMB spectrum allows for comparison with models 
that allow for departure from pure free-streaming (e.g., 
Cyr-Racine &: Sigurdson[[2012 ). The effect on the small- 
scaie power of a model with dark photons which are ini- 
tially coupled to dark matter and hence only start free- 
streaming after they decouple - implies that the phase 
shift and amplitude suppression associated with the free- 
streaming of radiation will not be uniform across all mul- 
tipoles. We leave the testing of such models to future 
work. 

4.3. Massive neutrinos 

In the previous subsection, we estimated the number 
of effective relativistic species, Aeff. If A^eS is assumed 
to arise solely from neutrinos impacting the CMB, these 
are also assumed to be massless in our fiducial model. 
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However, the CMB is sensitive to the sum of neutrino 
masses, which is related to the energy density of massive 
neutrinos via 



93.14 eV 



(11) 



In general there will be a degeneracy between the mass 
of a neutrino species and increased relativistic degrees of 
freedom at early times. After the neutrinos become non- 



WMAP7+ACT 

WMAP7+ACT+BA0+HST+ACTskewness 
WMAP7+ACT+BA0+HST 




0.4 0.6 

Zmv [eV] 

Fig. 8. — Marginalized one-dimensional likelihood for the sum of 
the neutrino mass, Sm^ in the case of the ACT+WMAP7 data, 
and with the addition of the BAO and Hq constraint. In addi- 
i :ion. we add a const raint on ag from ACT skewness measurements 
| |WiIson et al.|2012| l. 



relativistic, neutrino free streaming washes out structure 
on small scales. In Figure [8] we show the constraints on 
the sum of the neutrino masses assuming N^s is fixed to 
3.046. For the ACT in combination with WMAP7 data 
(and keeping the number of neutrinos fixed at N^s = 
3.046), we find 

Sm^ < 0.70 eV (WMAP7 + ACT, 95% CL). (12) 

Adding in distance measurements through the BAO and 
Hubble constant prior (which breaks the degeneracy be- 
tween YiTfiy and i?o)i improves the constraint to 

Em^ < 0.39 eV (WMAP7 + ACT + BAO+i7o, 95% CL). 

(13) 

Equations ((H]) and ([l3]) show that adding the ACT 
)]mr. 



data strengthens the const raints by factors of 1.5 and 1.8, 
respectively, compared to |Komatsu et al.^ ( 20111). These 
results are consistent with the |Hou et al ( |2012 1 upper 
limit of Yj-rrii, < 0.66 (95% CL) tor the ir combina t ion of 
CMB+BAO+i/o measurements. Whenl Hou et al.| ( |2012 | 
include SPT cluster measurements, this strengthens to 
a detection of non-zero neutrino mass, Etoj^ = 0.32 ± 
0.11 eV, with a corresponding value of as = 0.739 ± 
0.027; however, this value of as is in tension at the 3cr 
level with the AC T cluster mea surement of as = 0.829 ± 
0.024 JHasselfield et al | 120131 a nd the ACT skewness 
measurement ( Wilson et al.||2012| specified by Eq. 10 



Imposing the constraint on as from the ACT skewness 
measurement given in Eq. [TO] yields 

Sto^<0.30 eV 

(WMAP7 -I- ACT + ACTSkewness -t- BAO+Hq, 



95% CL). (14) 
Figure [9] illustrates the marginalized 68% and 95% con- 



WMAP7+ACT 

WMAP7+ACT+BA0+HST 

WMAP7+ACT+BA0+HST+ACTskewness 




Fig. 9. — Marginalized two-dimensional (68, 95%) contours in 
the N^ff — Sm^ plane. The number of effective relativistic de- 
grees of freedom is not strongly correlated with neutrino mass us- 
ing CMB data alone, however adding in complementary data sets 
helps tighten constraints. The star indicates the standard point 
(TVcfi, Em^) = (3.046,0.045 eV), where 0.045 eV is the minimum 
neutrino mass. 



tours for the sum of the neutrino masses and the Ngs 
when both are varied simultaneously. Adding in the Hub- 
ble constant and BAO data pushes N^s to slightly higher 
values, while the prior on as from ACT tSZ skewness 
measurements lowers the effective number of degrees of 
freedom. In all cases the sum of the neutrino masses is 
consistent with zero. 

4.4. Early dark energy 

The small-scale damping seen in the ACT data can also 
be interpreted as arising from through a non- negligible 
amount of dark energy at decoupling. The early dark 
energy (EDE) component may be specified through its 
density parameter (relative to the energy required for a 
fiat universe) Q^^ja) and an equation of state w{a), given 
by ( |Wetterich||1988| poran fc Robbers||2006[ ) 



w{a) = — 



r!dc(0) + r!™(0)a3-o 

1 (ilnf2(jc(a) 



3(1 — Odo(a)) rflna 3(a + Ccq) 



(15) 



The amounts of dark energy and matter today are 
given by r2dc(0) = flA,^m(fi) = flc + i^b respectively, 
and fleq is the scale factor at matter-radiation equality. 
D,c is the fraction of dark energy allowed at early times, 
while the present value of the equation of state of this 
early dark energy is expressed as wq- The scaling be- 
havior of the equation of state, tracking the dominant 
component at each cosmic era, gives rise to a radiation- 
like component at high redshifts, enhancing the damp- 
ing of the small-scale CMB tail. As the amount of EDE 
at early times tends to zero, this model approximates 
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TABLE 2 

ACDM AND EXTENDED 'DAMPING TAIL' PARAMETERS CONFIDENCE LIMITS FROM THE CMB DATA ALONE: ACT DATA 

COMBINED WITH SEVEN-YEAR WMAP DATA. 





Parameter"^ 


iior^iVi 




















-4- O 
-\- i fie 


1 /-V / f^n 


Primary 


^ 


2.251 ± 0.047 


2.237 ± 0.054 


2.242 ± 0.045 


2.257 ± 0.046 


2.260 ± 0.045 


A TV 

ACDM 




0.114 ± 0.005 


0.109 ± 0.011 


0.118 ± 0.006 


0.117 ± 0.006 


0.115 ± 0.005 




1009a 


1.040 ± 0.002 


1.041 ± 0.003 


1.040 ± 0.002 


1.059 ± 0.002 


1.060 ± 0.002 




ris 


U.y72 ± U.U12 


0.9d5 ± U.UzO 


r\ r\'7n i n r\-\ o 

0.970 ± U.012 


U.977 ± O.Uli 


0.97U ± U.013 




T 


U.Uyi it U.UlO 


u.uyi it U.UiO 


U.yio it U.Ulo 


U.UyU it U.UiO 


n QQQ -L n ni 
U.ooo it U.UiO 




l^fr/'i nlO A 2 \ 


o.iy ± U.U4 


Q on _L o n 
O.20 ± U.UO 


Q OA _l_ n r\n 
o.zO ± U.OO 


Q 1 o _L n n c; 
3.18 ± O.UO 


Q on _L n CiA 
o.zU ± U.U4 


Extended 






2.79 ± 0.56 










Em^ (eV) 






< 0.70 








Ce, UIQ 








< 0.03 






wo 








< -0.45 






a/ao 










1.004 ± 0.005 


Derived 




0.83 it 0.03 


0.81 ± 0.05 


0.77 ± 0.05 


0.72 ± 0.07 


0.84 ± 0.03 




Ha 


0.72 it 0.03 


0.72 ± 0.03 


0.69 ± 0.04 


0.62 ± 0.07 


0.73 ± 0.03 






0.28 it 0.03 


0.28 ± 0.03 


0.32 ± 0.05 


0.38 ± 0.07 


0.27 ± 0.03 




Ho (km s"l Mpc~^) 


70.0 it 2.4 


68.8 ± 3.5 


67.2 ± 3.1 


61.5 ± 5.0 


71.2 ± 3.6 


Secondary 


a-tsz 


3.4 it 1.4 


3.4 ± 1.4 


3.4 ± 1.3 


3.5 ± 1.4 


3.5 ± 1.4 




'ikSZ 


O.U 


^ o.^ 


^ o.o 


<^ o.^ 


^ S 1 




Op 


7.0 ±0.5 


7.0 ±0.5 


7.0 ± 0.5 


7.2 ± 0.5 


7.2 ±0.5 




Oc 


5.0 ±1.0 


5.0 ± 1.0 


5.0± 1.0 


5.2 ± 1.0 


5.1 ± 1.0 




as 


3.1 ±0.4 


3.1 ± 0.4 


3.1 ± 0.4 


3.1 ± 0.4 


3.1 ± 0.4 




/3c 


2.2 ± 0.1 


2.2 ±0.1 


2.2 ±0.1 


2.2 ± 0.1 


2.2 ±0.1 






0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 






0.4 ±0.2 


0.4 ±0.2 


0.4 ± 0.2 


0.4 ±0.2 


0.4 ± 0.2 


Calibration 


yls 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ±0.01 


1.02 ±0.01 


1.01 ± 0.01 




y2s 


1.03 ±0.02 


1.03 ±0.02 


1.03 ±0.02 


1.04 ±0.02 


1.04 ± 0.02 




yle 


1.00 ±0.01 


1.01 ± 0.01 


1.01 ±0.01 


1.02 ±0.01 


1.02 ±0.01 




y2e 


1.00 ±0.01 


1.00 ±0.01 


1.00 ±0.01 


1.00 ±0.02 


1.00 ±0.02 




-21n^ b 


8147 


8147 


8148 


8147 


8146 



^For one-tailed distributions, the upper 95% CL is given. For two-tailed distributions the 68% CL are showrn. 

'^The -2\nJff is based on the contribution from ACT (710 data points) and the WMAP7 UkeUhood (whose best-fit gives -2Inif=7478); 
including priors on the calibration parameters. The total degrees of freedom are therefore 8169 or 8168 if we consider the 19 or 20 parameter 
models. 

'^In Table 2 of |Dunkley et ar] | |2013| l, the secondary parameters, including calibrations, are reported assuming a best-fitting cosmological 
model. Here we estimate the calibration factors simultaneously with the ACDM parameters. 



TABLE 3 

Constraints on Ncb from the ACT data combined with seven-year WMAP data and other data. 



Data set 


A^cff 


Ho 




WMAP7+ACT 
WMAP7+ACT+SPT 
WMAP7+ACT+BAO 
WMAP7+ACT+HST 
WMAP7+ACT+BAO+HST 


2.79 ±0.56 
2.98 ±0.47 
2.87 ± 0.60 
3.18 ±0.46 
3.50 ±0.42 


68.8 ±3.5 
71.2 ±3.4 
67.7 ± 3.1 
72.2 ±2.2 
71.2 ±2.1 


0.109 ±0.011 
0.110 ±0.010 
0.114 ±0.011 
0.113 ±0.011 
0.124 ±0.009 



the standard ACDM cosmological model. We place a 
prior that w{a) > —1. This EDE model has been consid- 
ered by many authors (Ide Putter et a l.'2009'; HoUcnsteinl 
et al.||2009[ |de P utter et ai.| |!^010[ IL'alabr csc ct al. 2011; 

rese et al. 2011a b; Rcicha rdt et ar.||2012| showing 
however that current CMB observations combined with 
large scale structure data have no preference for a non- 
zero EDE d ensity. We include dark energy clustering as 
discussed in Calabrese et al. (2011) but we fix the dark 
energy sound speed and viscosity parameters to one and 
zero, respectively, as expected for a perfect fluid. Prior 
to this work, the most recent constraint on this model is 



from Reichardt et al. (2012) who report an upper limit 
of fie < 0.U18 at 95% conhdence from CMB only data, 
combining WMAP and SPT. In our analysis of the ACT 
data in combination with WMAP7 data and the ACT 
deflection measurement, we obtain the upper bound 

rie < 0.025 (WMAP7 + ACT + ACTDcfi, 95% CL). 

(16) 

with the bound on the equation of state found to be Wq < 
—0.45, see Figure [TOl We do not include the combination 
of early dark energy and Ncs , as the combination of these 
two parameters is largely unconstrained using the current 
small-scale CMB data. 
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TABLE 4 

ACDM AND EXTENDED SINGLE-MODEL INFLATIONARY AND DAMPING PARAMETERS FROM CMB DATA ALONE: ACT DATA 

COMBINED WITH SEVEN-YEAR WMAP DATA. 





Parameter^ 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 








LLlLg 1 yJL 111 t\j 


-t- r 


-1- Y 




-1- Oil 


Primary 


^ 

100f!(,/i 


2.251 ± 0.047 


2.253 ± 0.051 


2.248 ± 0.056 


2.236 ± 0.052 


2.246 ± 0.044 


2.279 ± 0.052 


ACDM 


L 2 


0.114 ± 0.005 


0.115 ± 0.006 


0.111 ± 0.006 


0.113 ± 0.005 


0.113 ± 0.047 


0.113 ± 0.005 






1.040 ± 0.002 


1.040 ± 0.002 


1.040 ± 0.002 


1.039 ± 0.003 


1.040 ± 0.002 


1.040 ± 0.002 




ris 


U.972 ± 


0.981 ± 0.031 


U.980 ± 0.U17 


r\ c\C'y 1 Ci r\-\ /i 

0.967 ± 0.014 


0.970 ± 0.011 


0.973 ± U.Olz 




T 


U.uyi it U.UiO 


U.uyz It U.Ui i 


U.uyz It U.UiO 


U.Uoy It U.Uio 


U.uyu It U.Uiu 


U.uyz it U.Uio 




i^.t/'i nlO A 2 \ 
log(10-^"A;^) 


o.iy ± U.U4 


3.19 ± O.Oo 


o.lD ± O.U/ 


0.2U ± O.Oo 


O.20 ± 0.04 


o.lo ± U.Uo 


Extended 


dns/d\n k 
r 

Yp 
G/i 




-0.004 ± 0.012 


< 0.30 


0.225 ± 0.034 


-0.031 ± 0.026 


< 1.7 X 10~' 


Derived 


0-8 


0.83 ± 0.03 


0.83 ± 0.03 


0.82 ± 0.03 


0.81 ± 0.03 


0.81 ± 0.03 


0.82 it 0.03 




Qa 


0.72 ± 0.03 


0.71 ± 0.03 


0.73 ± 0.03 


0.72 ± 0.03 


0.63 ± 0.08 


0.73 it 0.03 






0.28 ± 0.03 


0.29 ± 0.04 


0.27 ± 0.03 


0.28 ± 0.03 


0.40 ± 0.10 


0.27 it 0.03 




Ho (km Mpc~^) 


70.0 ± 2.4 


69.6 ± 2.9 


71.1 ± 3.1 


70.0 ± 2.3 


59.8 ± 8.0 


70.9 it 2.5 


Secondary 


atsz 


3.4 ± 1.4 


3.3 ± 1.4 


3.4 ± 1.4 


3.5 ± 1.4 


3.4 ± 1.4 


3.4 it 1.3 




<JkSZ 


s 

o.u 




^ o.o 




<^ O. 1 


u.y 




Op 


7.0 ±0.5 


7.1 ±0.5 


7.0 ± 0.5 


7.1 ±0.5 


7.0 ±0.5 


7.3 ±0.6 




file 


5.0 ± 1.0 


5.0 ±0.9 


5.0± 1.0 


5.0 ±1.0 


4.9 ±0.9 


4.9 ± 1.0 




as 


3.1 ±0.4 


3.0 ±0.5 


3.1 ± 0.4 


3.1 ±0.4 


3.1 ± 0.4 


3.1 ± 0.4 




/3c 


2.2 ±0.1 


2.2 ±0.1 


2.2 ±0.11 


2.2 ± 0.1 


2.2 ±0.1 


2.2 ±0.1 




Qge 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.2 






0.4 ± 0.2 


0.4 ± 0.2 


0.4 ±0.2 


0.4 ±0.2 


0.4 ±0.2 


0.4 ±0.2 


Calibration 


yls 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ±0.01 


1.02 ±0.01 




y2s 


1.03 ±0.02 


1.03 ±0.02 


1.03 ±0.02 


1.03 ±0.02 


1.03 ±0.02 


1.04 ±0.02 




yle 


1.00 ±0.01 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ± 0.01 


1.01 ± 0.01 




y2e 


1.00 ±0.01 


1.00 ±0.01 


1.00 ± 0.01 


1.00 ±0.02 


1.00 ±0.01 


1.00 ±0.02 




-2 In if b 


8147 


8147 


8148 


8147 


8148^^ 


8149 



^For one-tailed distributions, the upper 95% CL is given. For two-tailed distributions the 68% CL are shown. 

'^The -21nif is based on the contribution from ACT (710 data points) and the WMAP7 hkehhood (whose best-fit gives -2 lnif=7478); 
including priors on the calibration parameters. The total degrees of freedom are therefore 8169 or 8168 if we consider the 19 or 20 parameter 
models. 

'^In the case of the ACDM + model, both ACT power spectrum and deflection measurements are used, which increases the — 2i? value. 



WMAP7 

WMAP7+ACT 

WMAP7+ACT+Defl 




0.08 



Fig. 10. — One-dimensional marginalized likelihood for the early 
dark energy density, Oe. Adding in deflection measurements 
(shown by the inner-most curve) tightens the constraints on the 
allowed energy density in early dark energy. 



4.5. Inflationary parameters 

Inflation provides mechanism for the seeding of cosmo- 
logical structure through smaU fluctuations in the early 



universe. We constrain the spectral index of the initial 
spectrum of scalar density fluctuations through the mea- 
surement of the high— ^ tail of the angular power spec- 
trum. The scalar s pectrum of curvature pertu rbations is 
parameterized via ( Kosowsky &: Turner||1995 ) 



(17) 

We constrain the amplitude A|j(A:o), spectral index 
ns{ko) and 'running' of the spectrum dus/dlnk, defined 
at a pivot point k = 0.002Mpc~^. The amplitude of 
scalar perturbations is found to be 

log(10i°A^) = 3.19 ± 0.04 (WMAP7 + ACT), (18) 

where the factor of 10^" is a numerical factor included 
to ensure robustness to numerical precision errors while 
sampling. 

4.5.1. The scalar spectral index 

A generic prediction of inflationary models is a nearly 
scale invariant spectrum; any de yiations provide power- 
ful te sts of inflationary models fMukhan ov fc Chibisov 
198T|[Hawking.l982,,Starobinsky.l982,,Gutli P i| 
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TABLE 5 

Extended model parameters confidence limits from the ACT data combined with seven-year WMAP data, 

BAO and Ho measurements. 





Parameter^ 



















1 AT 


+mi, 


V 




Primary 


lOO^fc/i^ 


2.225 ± 0.038 


2.271 it 0.046 


2.259 it 0.039 


2.256 ±0.040 


2.261 it 0.044 


ACDM 




0.116 ± 0.003 


0.125 it 0.008 


0.115 it 0.002 


0.115 ±0.003 


0.114 it 0.005 




1009a 


1.040 ± 0.002 


1.038 it 0.002 


1.040 it 0.002 


1.040 ± 0.002 


1.040 it 0.002 




Us 


0.971 ± 0.009 


0.982 ± 0.014 


0.974 ± 0.010 


0.972 ± 0.010 


0.974 ± 0.011 




T 


0.090 it 0.014 


0.092 it 0.015 


0.091 it 0.014 


0.089 it 0.014 


0.092 it 0.015 




Iog(10iOA|,) 


3.20 it 0.03 


3.19 it 0.03 


3.19 it 0.03 


3.19 it 0.04 


3.19 it 0.04 


Extended 






3.50 it 0.42 










Sm„ (eV) 






< 0.40 








r 








< 0.16 
















-0.0020 it 0.0047 


Derived 




0.83 it 0.02 


0.86 it 0.03 


0.79 it 0.04 


0.83 it 0.02 


0.83 it 0.03 




Ha 


0.71 it 0.01 


0.71 it 0.01 


0.71 it 0.01 


0.71 it 0.01 


0.72 it 0.01 






0.29 it 0.01 


0.29 it 0.01 


0.29 it 0.01 


0.29 it 0.01 


0.29 it 0.01 




Ho (km s-l Mpc-i) 


69.3 it 1.0 


71.2 it 2.0 


68.8 it 1.0 


69.3 it 1.0 


69.0 it 1.1 


Secondary 


ItSZ 


3.4 it 1.4 


3.4it 1.4 


3.4 it 1.4 


3.4 it 1.4 


3.4 it 1.5 




"kSZ 


< 8.3 


< 8.8 


< 8.5 


< 8.2 


< 8.5 




dp 


7.0 it 0.5 


7.0 it 0.5 


7.0 it 0.5 


7.1 it 0.5 


7.0 it 0.6 




Oc 


5.0 it 1.0 


4.8 it 1.0 


5.0 it 1.0 


5.0 it 1.0 


5.0 it 1.0 




Os 


3.1 it 0.4 


3.0 it 0.4 


3.1 it 0.4 


3.1 it 0.4 


3.1 it 0.4 






2.2 it 0.1 


2.2 it 0.1 


2.2 it 0.1 


2.2 it 0.1 


2.2 it 0.1 




Oge 


0.9 it 0.2 


0.9 it 0.2 


0.9 it 0.2 


0.9 it 0.2 


0.9 it 0.2 






0.4 it 0.2 


0.4 it 0.2 


0.4 it 0.2 


0.4 it 0.2 


0.4 it 0.2 


Calibration 


yls 


1.01 it 0.01 


1.01 it 0.01 


1.01 it 0.01 


1.01 it 0.01 


1.01 it 0.01 




y2s 


1.03 it 0.02 


1.03 it 0.02 


1.03 it 0.02 


1.03 it 0.02 


1.03 it 0.02 




yle 


1.01 it 0.01 


1.01 it 0.01 


1.01 it 0.01 


1.01 it 0.01 


1.00 it 0.01 




y2e 


1.00 it 0.02 


1.00 it 0.01 


1.00 it 0.01 


1.00 it 0.01 


1.00 it 0.01 




-21nif b 


8153 


8152 


8153 


8153 


8152 



^For one-tailed distributions, the upper 95% CL is given. For two-tailed distributions the 68% CL are shown. 

''The -2\nJC is based on the contribution from ACT (710 data points) and the WMAP7 hkehhood (whose best-fit gives -2 lnif=7478); 
including priors on the calibration parameters. The total degrees of freedom are therefore 8169 or 8168 if we consider the 19 or 20 parameter 
models. 



Bardeen, Steinhardt, fc Tumer|[l983l |Mukhanov, Feld-| 
man, & Brandenberger"1992|. Deviations from scale 



invariance have been te sted m a variety of mod els and 



contexts (e.g., Wang e t al.||199 9; Tcgmark & Zaldar 
riagal2002tlBrrdlc ct al. MI^THanncstad 2003, .Martin fc 



ningeval 2004a b; Scalfon ct al. 2005; Spe rgel et al.|2007[ 
Verde & Peirisi|2008| i Peiris fc Verdei2010nVazquez et al'| 



Peiris||2008l |Pe 
nth the hign-re 



^ ^ Lepz 

2011 1, with the higt-resolution of ACT expanding pos- 
sibilities for direct measurements of the power spectrum 
( Hlozek et al.|2012| ). The standard cosmological scenario 
in which no variation of the spectral index is allowed (i.e. 
dns/d\nk = 0) provides an excellent fit to the data. We 
obtain a constraint on the scalar spectral index of 



0.972 ± 0.012 (WMAP7 + ACT). 



(19) 



When dus/dlnk is allowed to float, moves to slightly 
different values as shown in Figure [TT] When considering 
the ACT-S and ACT-E spectra separately (while still in 
combination with WMAP7), wc find 



ns(ACT- E) = 
n,(ACT-S) = 



= 0.981 ±0.012; 
= 0.966 ±0.012. 



(20) 



The ACT-E data specifically prefer a higher value of 
the ACT-S data, wh ich is consistent with U g = 0.962 ± 
0.013, the value in Dunkley et al. (20111. However, 



the cosmological models that describe the two statist! 



cally independent data sets are consistent as shown in 
Appendix IXj The difference in the marginalized Ug is 
indicative of residual correlations between parameters. 
Combining the two data sets increases the value of Ug 
relative to WMAP7 alone. It is similarly higher than 
rig = 0.9538 ± 0.0081 from WMAP7+SPT+BAO+H0 



reported in Story et al. (2012 1. 
Including BAU and Hq data yields the constraint 

ug = 0.971 ± 0.009 (WMAP7 + ACT + BAO + Ho), 

(21) 

which rules out a scale-invariant Ug = 1 spectrum at a 
significance of 3.2cr. The ACT+WMAP7 value of Ug = 
0.972 ± 0.012 is ~ 1.2a higher than t he recently released 
WMAP9 result n, = 0.9579±0.0081 (IBennett et al.|2012 



Hinshaw et al. 
SP'i' data and 



2012), however this constraint includes 



he improved WMAP data. 



4.5.2. ACDM + Running Index 

In addition to the spectral index alone, we test for de- 
viations from a perfect power law, in the form of running 
of the spectral index. Our cosmological models are run 



around a pivot point of ko = 0.002 Mpc . Figure 11 
shows the joint constraints on the spectral index and its 
nmning at a de-correlated pivot point fco — 0.015 Mpc^^, 
chosen to minimize the corre lation between the two pa- 
rameters (Cortes et al. 2007). We show the constraints 
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1.00 



^ 0.98 



0.96 




WMAP7 
WMAP7+ACT 



1.10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 

diis/d \nk 

Fig. 11. — Two-dimensional marginalized limits (68% and 95%) 
for the spectral index, ris, plotted at the pivot point fcg = 
0.015 Mpc~^, and the running of the index dus/dlnk, for 
ACT+WMAP7, compared to WMAP7 data alone. The data are 
consistent with the model dus/dlnk = 0. 



for the ACT data combination with WMAP7, compared 
to the WMAP7 data alone. This relation between the 
indices at these two pivot points is given by: 

ns{ko = 0.015 Mpc"^) = r7.s(fco = 0.002 Mpc"^) 

+ ln(0.015/0.002)- ' 



din fc 



(22) 

The ACT data are consistent with no running of the 
spectral index, with 



dug 
din fc 



-0.004 ± 0.012 (WMAP7 -I- ACT). (23) 



4.5.3. A CD M+ Tensor modes 

The standard ACDM model assumes that the initial 
fluctuations are purely scalar modes; however tensor fluc- 
tuations can also be seeded at early times, generating 
gravitational waves. While gravitational waves have not 
yet been detected, limits can be placed on the power of 
tensor fluctuations by constraining the so-called tensor- 
to-scalar ratio 



r = A)l(fco)/A|(A;o), 



(24) 



where A^(fco) is the amplitude of primordial gravi- 
tational waves, measured at a reference scale fco = 
0.002 Mpc^^; these tensor fluctuations contribute to the 
CMB in both temperature and polarization. Small-scale 
temperature measurements such as those from ACT help 
to constrain the amplitude r by reducing its degeneracy 
with the scalar spectral index, Ug. 

This degeneracy arises as models with large values of 
r have more power from tensor modes at small £ and 
therefore require smaller values of the spectral index to 
compensate. The analysis of the 300 square degree ACT 
1-year data combined with WMAP7 data constrained 
r < 0.25 (95% CL) using only CMB data. In this 3-year 
analysis we find the limit on r has increased, with 





WMAP7+ACT 

WMAP7 + ACT + BAO + HST 



Fig. 12. — Top pane/.Two-dimensional marginalized distribution 
(68% and 95% CL) for the tensor-to-scalar ration r and the scalar 
spectral index ris. The brown contours show the constraints using 
only WMAP7 data; the blue contours show the constraints when 
adding the 3-year ACT data, and BAO and Hq measurements, 
and constrain r < 0.16. Recfast vl.5 was used to compute the 
constraints. Bottom panel: One-dimensional marginalized con- 
tours for the tensor-to-scalar ratio for the ACT data in combi- 
nation with WMAP7 (dashed curve) and for the combination of 
WMAP7+ACT+BAO-l-iTo (solid curve). 



We find an upper bound w 1 5% higher for the ten sor-to- 
scalar ratio than reported in Dunkley et al. ( 2011[ ). This 
is driven by the Ug — r degeiieracy and the tact that the 
ACT 3-year constraint on Ug is higher than the previ- 
ously reported value. The curvature of the mar gina lized 
likelihood is shown in the bottom panel of Figure 12 when 
only using CMB data. Once distance measurements are 
added the bound tightens. It is worth noting that the 
WMAP7-only constraints on the tensor-to-scalar ratio 
and sca lar spectral index als o differ from the published 
results jK omatsu et al.|20TT ) , as can be see n by compar 



ing Fi gure Tl in this work to Figure 13 in Larson et al 



(2011 1. These arise from changes in the version of Rec 



r < 0.30 (WMAP7 + ACT, 95% CL). 



(25) 



fast used in the WMAP7 version of this particular pa- 
rameter combination. These differences are summarized 
in Table m 

When adding in data from BAO and the Hubble con- 
stant value, we obtain 

r < 0.16 (WMAP7 + ACT + BAO+i7o, 95% CL). (26) 

This value is generally consistent with the recently re- 
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TABLE 6 

Constraints on the spectral index and tensor-to-scalar ratio for different versions of the Recfast code. 



Data set 


Recfast version 


ris 


r (95% CL). 


WMAP7 alone 


vl.4.2 


0.981 ± 0.020 


< 0.32 




vl.5 


0.990 ±0.021 


< 0.39 


WMAP7 -1- ACT 


vl.5 


0.980 ±0.017 


< 0.30 



leased WM AP9 results ( [Bennett et al.||2012[ Hinshaw 
et al.||2012[ ), which find a tighter bound using WMAF9 
data m combination with smaU-scale CMB spectra, BAO 
measurements, and the Hubble constant measurement, of 
r < 0.13. 



4.6. Curvature 

One of the simplest extensions of the primary cosmo- 
logical model is allowing the curvature of the universe to 
vary, through the relation 



(27) 



where fl^ is the value of the curvature at z = 0. A neg- 
ative sign for fl^ implies a c losed universe, while flh > 



describes an open universe. Guth & Nomura (2012) de 



scribe how a measurement of non-zero curvature places 
limits on not only the geometry of the universe, but also 
on the space of allowed models of inflation; should pos- 
itive curvature be detected at the level of fl^ < —10^^, 
the current paradigm of eternal inflation would be falsi- 
fied, as this scenario generically predicts positive values 
of Ofe. A non-zero measurement of negative curvature 
(Ofc > 0) allows one instead to place constraints on the 
different possible pre-inflationary histories and the prob- 
ability measures within those frameworks. 

When combining the CMB data with measurements of 
the Hubble constant and BAO we constrain the curvature 
to be 

Ofc = -0.0020 ±0.0047 (WMAP7 + ACT + BAO+i?o), 

(28) 

a value consistent with a flat universe, and the previ- 
ously reported WMA P7 resuh of flk = -0.0023;5j;EJ^^6 
( [Komatsu et al.||2011| . 



4.7. ACDM + the fine structure constant a 

Given the damping of the power spectrum observed 
at large multipoles, one might ask if this damping can 
be explained in the context of models which allow for 
the variation of fundamental constants, such as the fine 
structure constant, a = e'^/h. Dark energy is a possible 
source of variation of the fine structure constant, and 
hence limits on the temporal variation i n a put limits on 
the e quation of state at late times (e.g., Parkinson et al. 
20041. The parameter a mediates the strength of the 
electromagnetic interaction, and therefore it affects the 
formation of CMB anisotropies by changing the scatter- 
ing rate of the Thomson scattering process as 

f = XeUeCaT, (29) 



: 87r(Q;/me)^/3 is the Thomson cross section 
c ~ 1 and Xe is the electron fraction, itself depen- 



where ctt 
for h 

dent on the scale factor of the universe, a{t) (e.g., Battye 



et al.|[200T and references therein). 




1.00 1.02 

a/oo 

Fig. 13. — One-dimensional marginalized likelihood for a/ao 
from the ACT data in combination with WMAP7. The constraints 
arc consistent with no variation in a at recombination. 



We quantify this physical mechanism in terms of the 
ratio a/ao, where an = 1 /(137.036) is the standard local 
value (Mohr et al. 2012) and a is the value during the 



recombination process. A value a/ao > 1 will generate a 
shift in the acoustic peaks towards higher multipoles as 
well as less damping. M cncgoni ct al. (2012) allowed for 
a variation in a and illustrated the degeneracy between 
such variations and, for example, A^effj at early times. 
If variations in a arise because of some coupling of the 
photon to a scalar field ( Bekenstein|2002 1 , in general one 
might e xpect fluctuations in a to be both temporal and 
spatial ("Sigurd son et al.|2 003) , however for simplicity we 
ignore spatial variations here. In Figure [TSl we show the 
constraints on the variation in a for the^CT+WMAP7 
data combination, which give a/ao — 1.004± 0.005, con- 
sistent with the standard picture. 

4.8. Lensing 

The lensing of the CMB photons by gravitational struc- 
ture smears out the peaks and troughs of the CMB tem- 
perature spectrum, transferring power between modes. 
AC T reported the fi rst detection of intrinsic CMB lens- 



ing ( Das et ar[2011a ) at 4cr; in many models this provides 
evidence tor dark energy from the CMB alone (Sherwin 



et al.|20lT ) , adding further support to the ACDM model, 
i'he youth Pole Telescope (SPT) presented a det ection of 
the lensing defiec tion power spectrum at 6.3(t (van En- 
gelen et al. 2012). The strength of the lensing can be 
simply parameterized through a lensing amplitude 
( [Calabrese et al.||2008[ ) 

(30) 



ArCt 



theory 



where Cf 
and Al = 



is the power spectrum of the lensing field 0, 
1 corresponds to the amount of lensing in the 
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1.0 



0.8 



§0.6 



03 



0.2 




WMAP7+ACT 
WMAP7+ACT+Defl 



2.5 



3.0 



3.5 



relation on the helium abundance. Here we also con- 
sider the case where both parameters are allowed to vary 
freely. The constraints fro m t his model are shown in the 
right-hand panel of Figure [15] with the BBN consistency 
relation shown as a line in the plane. The model with 
Ngff = 3.046 and Yp = 0.24 is shown as a star. 

4.10. Cosmic strings 

While cosmic strings are no longer a viable seed of 
all cosmic structure, should a sub-dominant connected 
gas of cosmic strings exist, they would add power to the 
small-scale tail of the CMB power spectrum. The am- 
plitude of diffuse power is proportional to the string ten- 
sion, G^, and hence small-scale measurements provide a 
handle to constrain various string models. While direct 
searches for map linear string signatures in CMB maps 



Fig. 14. — One-dimensional marginalized likelihood for the lens- 
ing parameter from WMAP7 and ACT temperature and deflection 
spectrum measurements. 



standard ACDM model. As this parameter multiplies the 
entire lensing power spectrum, it is defined over the same 
range in multipole space used for the theoretical temper- 
ature spectrum. We measure the lensing amplitude from 
the smearing of the temperature power spectrum as 



provide interesting constraints on the string tension (Be- 
vis et al.||2010 Mukherjee et al.|["2011 ), the simplest pre- 



scription toFmcIudTngTEie~diIIiise power is to model it as 
a network of string segments which produce loops that 
decay into radiation. 



We consider a N ambu string template as used in Bat 



tye & Moss ( 2010 1, modeled as a power law invers ely pro 



portionaTTo multipole on especially small scales (Fraisse 



et al. 2008). Various authors mo del the st ring net 



work iismg different approaches (e. g., Bennett &i Boucl^t] 



1.7 ±0.38 (WMAP7 + ACT). (31) [l990| [Pogosian fc VachaspatHl999[ [Battye fc Moss|20 



When adding in the deflection spectrum measurement, 
this value is reduced to 

Al = 1.34 ± 0.23 (WMAP7 + ACT + ACDefi), (32) 

as is shown in Figure [T4j While this value is somewhat 
high, it is consistent with A^, = 1 at the 1.5a level, where 
Al = 1 is the amplitude of lensing in the ACDM sce- 
nario. 

4.9. Primordial Helium 

The conditions at decoupling are highly sensitive to 
the expansion rate of the universe at early times. Ex- 
tra relativistic species change the expansion rate of the 
universe at a given temperature. An increased expan- 
sion rate cools the universe faster, and nuclear freeze- 
out occurs at an earlier time (recall that H = l/{2t) 
at early times), when the concentration of neutrons is 
higher, changing the primordial abundances of light el- 
ements including ^He. We parameterize the primordial 
helium mass fraction as 

Yp = 2x, (33) 

where x = nn/{nn + rip) is the neutron abundance, n„ 
is the number densi ty of neutron s and rip the number 



where the equations of motion are solved using either 
Nambu or Abelian-Hig gs methods. Other app roaches al- 



low for more freedom ( [Foreman et al. 2011), including 
the mean velocity of strings, correlation length of the 
string network, and small-scale structure of the strings 
as parameters in the model. 

Previous constraints on the string tension from the 
ACT 1-year data in combination with WM AP7 limited 



G^l < 1.6 X 10"^ (Dunkley et al. 2011), while con- 



straints from the WMAP7 data m combination with SPT 



( Dvorkin et al.| 2011) using a "Nambu-Goto" method 
bound the allowed string tension to be G/i < 1 x 
10^^ (95% CL). In this analysis, we constrain 

Gfi < 1.15 X lO"'^ (CMB + Dcfl + BAG + Ho, 95% CL). 

(35) 

Early analyses of CMB data at {i < 1000) in the 
context of string models found that the model with a 
Harrison-Zeld'ovich spectrum with = 1 and cosmic 



strings was consiste nt with the data (Bevis et al. 2008 



Battye et al. 2006) (and even that the six-parameter 



model with both strings and a scale-invariant pertur 
bation spectrum fit the data as well as the ACDM 
model without strings), while more recent studies have 
shown that this model no longer p rovides an ade- 



density of protons (Peebles 1993). The abundance of 
primordial helium is degenerate with Ncs in a model with 
standard Big Bang Nucleosynthesis through Eq. ([8[). The 
constraints from WMAP7 in combination with ACT data 



are: 



Yp = 0.225 ± 0.034 (WMAP7 + ACT), (34) 

where N^ff is held fixed at the standard value of N^g = 
3.046. The marginalized one-dimensional likelihood in 
this case is shown in the left-hand panel of Figure [TS] 



When constraining N^s in general, as in Section 4.2, 
and Tables |2] and [^ we impose the BBN consistency 



quate fit to small-scale CMB data (Urrestilla et al. 
2011[ ). We find that fixing the primordial power 
to be scale invariant, and combining all CMB data 
(WMAP7-hACT-hSPT-HACTDefl) with BAG and the 
Hq measurement yields a constraint on the string am- 
plitude of qstr = 0.005, or G^i = 7.1 x 10"^ (95% CL). 
The difference in effective values between the two 
models [rig = 1 + g^tr versus rig free and no strings; 
both models have the same number of parameters) is 
-21n^ACDM + 21n.if„^=i_strings = 32, indicating an 
5.7(7 preference for the ACDM model over a model with 
strings (using a Nambu string template) and a scale- 
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Fig. 15. — Left panel: Marginalized one-dimensional constraints on the helium fraction Yp while fixing the number of effective relativistic 
degrees of freedom, Ngff = 3.046. The dashed line indicates the standard value of = 0.24. Right panel: Marginalized two-dimensional 
68% and 95% contours for the effective number of relativistic degrees of freedom and the primordial helium abundance. The line indicates 
the BBN consistency relation (given in Eq. jsjl), while the star labels the ACDM model with N^s = 3.046 and Yp = 0.24. 



invaria nt spectrum. Recent work by [Lizarraga at al. | 
(2012) note that the string amphtude is highly correlated 
with iVoff, with A'cff being pushed towards much higher 
values when both the string tension and the number of 
effective degrees of freedom are varied together. 

5. SECONDARY PARAMETERS 

The assumed priors place d on the secondary p arame- 
ters are discussed in detail in Dunkley et al.| ( 2013[ ), where 
the values for assumed priors are presented and justified. 
We summarize and interpret our constraints here. 

5.1. Point source power 

The power spectrum includes a contribution from point 
sources that lie below our detection threshold. At 148 
GHz and 218 GHz, both synchrotron emission from radio 
galaxies and infrared emission from dusty galaxies con- 
tribute to the observed power. Point sources detected in 
either band are masked down to a flux limit of 15 mJy. 
We impose a prior of = 2.9 ± 0.4 at I = 3000 on 
the residual power from radio galaxies at 148 GH z from 
sources below this threshold (Gr alla et al. 2013), based 
on models fror n . Tucci et al. (2U11). I'br comparison, Re- 
ichar dt et al. (2012) use a model from De Zotti et al. 
( |20U5p , which would predict a residual power below 15 
mJy that differs b y about 5% from the prediction of the 
Tucci et al. (2011 ) model, and so lies within the Icr limit 
of the prior we impose. 

For the ACDM model, the error on the unresolved ra- 
dio galaxy power (at 150 GHz) is dominated by the prior 
on the amplitude of power a,t £ = 3000, with 



3.1 ±0.4. 



(36) 



We constrain the clustered (ac) and Poisson (op) power 
from GIB sources respectively to be 



5.0 ± 1.0; 



7.0 ±0.5, 



(37) 



where again, the amplitudes are for templates normalized 
to 1 /LtK2 at £ = 3000 and 150 GHz. 

The frequency dependence of the GIB power is modeled 
as the square of a modified blackbody (see Section 2.3 



of Dunkley et al. 20131, with effective dust temperature 



fixed to 1,1 = 9.7 K. The effective emissivity index is 
constrained to be /3c = 2.2 ±0.1. The corresponding 
constraint on the GIB spectral index between 150 and 
220 GHz is ad = 3.68 ± 0.10, where ad is defined in 
flux density units as 5*^ cx i'"''. This result is in good 
agreement with the earlier AGT 1 -year constraint of ad — 
3.69 ± 0.14 (Dunkley et al. |2011|), and somewhat higher 
than, though not inconsistent with, ad — 3.56 ± 0.07 
(or 3.45 ±0.11, w hen a tSZ-CIB correla tion is allowed), 
reported by SPT ( [Reichardt et al.||2012[ ). 

While /3c cannot be directly associated with the emis- 
sivity of individual GIB sources, we expect the GIB 
anisotropy measured by AGT to improve dusty star- 
forming galaxy constraints when used in conjunction 
with d ata at higher frequencies (as in Addison et al. 
2012a), and through ongoing cross-correlation analyses. 



5.2. Thermal SZ 

The amplitude of the thermal SZ (tSZ) power spec 
trum is a strong function of the ampl itude of matter 
fluctuations, and is very sensitiv e to dKomatsu & Ki 

8.3 

8 



tayama" 19991) 
et al. 2011|~_ 



scaling Hke o-f -^ (^Shaw et al. 2010 Trac 



We mode l the tSZ power using a template 



(Battaglia et al. 2011) from hydrodynamic simulations 
of cosmological volumes that includes radiative cooling, 
star formation, and AGN feedback. For erg — 0.8 this 
template has a predicted amplitude of Otsz = 5.6 at 
t — 3000. In this analysis of the AGT 3-year d ata we 
separate the t hermal and kinetic SZ components (Dunk- 
ley et al. 2013 1 into two dimensionless spectra normalized 
to unity at £o = 3000 in Eq. |3] Hence, the theoretical 
expectation for the tSZ amplitude at this £o = 3000 is 
5.6. We constrain the tSZ amplitude to be 



atsz =3.4±1.4. 



(38) 



Our current result corresponds to a measurement of erg — 
0.75lg'Q5 and is consistent with our prev iously reported 



values from atsz of cg = 0.77 ± 0.05 in [Dunkley et al. 
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( 2011 1 (which assumed a CTg scahng), and from ACT mea- 
surements of the th ree point function, erg = 0.79 ± 0.03 
(Wilson et al. 2012). Our tSZ constraint is also consis- 
tent with the SPT-rep orted value of Otsz = 3.3 ± 1.1 
( [Reichardt et al.||2012^ 



In the AC'l' a nalysis using only the 1-year data ( Dunk- 
ley et al. 2011 1 we modeled the SZ power spectrum us- 
ing only a total (dimensionless) amplitude of the SZ 
signal (thermal plus kinetic), which was constrained 
relative to a theoretical spectrum with erg = 0.8 as 
asz/asz,<Ts=o.8 = 0.85 ± 0.26 using both 148 GHz and 
218 GHz frequency bands, and which can be interpreted 
as i(e+l)/2TTCl°^^^ = 6. 8 ± 2.9 /^K2 at ^ = 3000 using 



the Battaglia et al. (2010) template. 

However, the connection between a tSZ amplitude and 
constraints on (jg depends on the modeling of the clus- 
ter physics. For a fixed value of erg the predicted am- 
plitude of the tSZ power spectrum at £ = 3000 varies 
by tens of percent depending on the tSZ model (e.g . , 



Battaglia et al.|2010[plaw et al.|2010||Trac et al.|2011[ ). 

All the direct measurements of cluster plwsics are limited 



to high-mass, low- redshift clusters (e.g. Planck CoUabo- 
ration et al. 2012), which do not constrain tSZ models, 
since around halt of the contribution to tSZ amplitude 
at £ = 3000 comes from lower mass and higher redshift 
cluste rs and groups (e.g., Trac et al. 2011 Battagha et al. 



The uncertainty in the tSZ model leads to a system- 
atic error in erg, which is not included in the cg error 
bar. Combining the power spectrum measurements with 
higher order correlations and cluster number counts will 
provide tighter constraints on erg and the tSZ models 
( [Bhattacharya et al.||2012[ [Hill fc Sherwin|[2012l ). 



5.3. tSZ-CIB correlations 

The tSZ-ClB correlation term is important in inter- 
preting the power from the kSZ effect, as the frequency 
dependence of the tSZ-CIB and kSZ are similar across 
the ACT and SPT freque ncy channels (see e.g., Figure 
5 of Addison et al.|[2012b ) . A larger value of ^ removes 
more power at H. ~ 3000, which can be compensated by 
an increase in the kSZ amplitude. As shown in Figure [3] 
^ is largely unconstrained by the ACT data, and the 
choice of the ^ < 0.2 prior has little impact on the cos- 
mological analysis. Adopting a wider prior of ^ < 0.5, 
estimated as the upper limit from the modeling by |Addi^ 
sonet al. (2012b I, degrades the kSZ amplitude constraint 
trom Oksz < 8.6 to aksz < 9.4 (95% CL), shown by the 
unfilled, dashed contours in Figure [16] The kSZ con- 
straint from fixing ^ = is Oks z < 8.0, indicated by the 
dotted horizontal line in Figure 16 The anti-correlation 
between the tSZ-ClB and kSZ power is more apparent in 
the joint fit with the SPT spectra. 

5.4. Kinematic SZ 

The kSZ power spectrum can be separated into two 
components. The first is a low-redshift component 
(which we refer to as the homogenous kSZ) that is 
sourced by perturbations in the free electron density after 
reionization with peculiar velocities; on large s cales this 
is known as the Ostriker-Vishniac (OV) effect ([Ostriker 
FVShmac]|l986t 

1970p . Tlie second is a high-redshift component that is 



WMAP7+ACT+SPT, ^ < 0.5 
WMAP7+ACT, 5 < 0.5 
WMAP7+ACT, ^ < 0.2 , ' ' 




-4 -3 -2 

Fig. 16. — Marginalized two-dimensional 68% and 95% con- 
tours for the amplitude of the kSZ power and the power in the 
tSZ-CIB correlation, for the combinations of ACT+WMAP7 and 
ACT-I-WMAP7+SPT. The filled contours are for the prior ^ < 0.2, 
while the unfilled contours show the degrading of constraints when 
weakening the prior to ^ < 0.5. The dotted line shows the limit on 
ikSZ for € = 0. 
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ACT 

ACT+SPT 

ACT+SPT,|<0.5" 

SPT 

SPT, I < 0.5 




Fig. 17. — Constraints in the z — plane from the patchy 
reionization signal at £ = 3000. The solid curves are for the ACT 
data in combination with WMAP7 data, the light (dark) dotted 
green curves are for SPT-fWMAP7, under the assumption of priors 
on the tSZ-CIB correlation of g < 0.5 < 0.2). The light (dark) 
dashed blue lines show the combination for WMAP7 data, ACT 
and SPT, also under the priors of § < 0.5 (^ < 0.2); the allowed 
regions are below t he colored lines. The g rav region is excluded 
by Gunn- Peterson ( |Gunn fc Peterson||1965 I observations and the 
WMAP7 limits on z are shown by the dark gray (Icr) and black 
(2(t) lines. 



sourced by fluctuations in the ionized fraction as well as 
the electron de nsity during reionization, which w e refer to 
as patchy kSZ ( |Knox et al.]|1998| [Gruzinov fc Hu 1998). 



as a function of £ (e.p 


■ 1 


McQuinn et al. 112007 


Zahn et al. 


2012 Mesinger et al. 

' rni ' ' 1 ' ' 


2012). 



tion with WMAP7 data is 

akSz < 8.6 (WMAP7 + ACT, C < 0.2, 95% CL). 



(39) 



When relaxing the prior on the tSZ-CIB correlation 
to < ^ < 0.5 and including SPT data, we obtain the 
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bound 



aksz < 6.9 (WMAP7 + ACT 
C < 0.5,95% CL). 



SPT, 



(40) 



We approximate the homogeneous kSZ contribution at 
i — 3000, aks z.hom = 1.5 ± . 0.5 using the scaling rela- 
tion from Shaw et al. (2012 1, however we do not include 
the uncerlainty m our analysis. We remove this homo- 
geneous kSZ contribution leaving only the contribution 
from patchy kSZ, aksz, patchy < 7.1. The amplitude of 
the patchy kSZ places constraints on the mean redshift 
(z) and duration of (A^), reio nization via the relation 
given by IBattaglia et al.l (|2012|) as 



OkSZ, patchy 



2.02 



1 



11 



0.12 



1.05 



0.47 



(41) 



The curves in Figure [17] illustrate the constraints on 
z and A, derived from combinations of ACT, SPT, and 



WMAP data assuming the Battaglia et al. (2012) scal- 
ing relation. The grey and black lines show the 68% and 
95% confidence regions from large -scale WMAP 7- yea r 



CMB polarization measurements (Larson et al. 2011). 



The ACT data mildly prefer a higher kSZ amplitude than 
SPT, which leads to an increase in the kSZ upp er li mit for 
WMAP7-t-ACT+SPT (dashed fine in Figure [l7|, com- 
pared to the WMAP7-I-SPT case (dotted fine). The 
choice of ^ prior drives the constraint when SPT data 
are included; constraints on the tSZ-CIB correlation and 
thus kSZ power should be significantly improved with 
cross-correlations of the ACT maps with CIB-dominated 



Herschel/SPIRE maps (Addison et aL][2012b |. 

Recently |Zahn et al.| ( 12012] ) obtained < 7.9 at 95% 
confidence iasing SP'l' data, although their analysis used 
a somewhat different treatment of the tSZ-CIB corre- 
lation, definition of A^ and reionization model to that 
considered here. 



6. CONCLUSIONS 

We have presented parameter constraints from the 
ACT 3-year data in the ACT-S and ACT-E regions, over 
three seasons, at 148 GHz and 218 GHz. Using the multi - 
frequency hkelihood presented in Dunkley et al. (2013), 
the data provide constraints on both the primary cosmo- 
logical model and a suite of secondary parameters. We 
confirm that the ACT data alone are fully consistent with 
the standard ACDM cosmological model. Our cosmolog- 
ical model is also consistent with data from the South 



Pole Telescope ( ]Reichardt et al.|20l"2||Keisler et al.|2011] ) 

While this work was being completed, updated c osmolog- 



al. 



ibina 



ical results were presented b y the WMAP team (]Bennett 
et al.]]2012| [Hinshaw et" 
the AC'l' data in comt 



2012) . Our results, based on 



ion with WMAP7 data, are 
consistent with the WMAP9 results, which have been 
tightened through the addition of two years of additional 



data and improved analysis techniques. We discuss gen- 
eral consistency between data sets in Appendix [X] 

We measure a power spectrum from the thermal 
Sunyaev-Zel'dovich effect that is consistent with ACDM, 
and place an upper bound on the power from the kine- 
matic SZ effect. This amplitude provides constraints 
on the time and duration of reionization; future cross- 
correlation studies will improve this bound. 

When combining the ACT data with data from the 
seven-year release from the WMAP satellite, we con- 
strain a host of extended cosmological models. We see no 
evidence for mean spatial curvature, running of the scalar 
perturbation spectral index, or extra relativistic degrees 
of freedom beyond three light neutrinos. The amount of 
primordial ^He is also constrained by the ACT data, and 
is consistent with the standard prediction from Big Bang 
Nucleosynthesis. 

The ACT data tighten the bound on the sum of light 
neutrino masses to 0.39 eV, using a measurement of trg 
from the skewness of the ACT temperature map. We 
bound the allowed contribution of tensor modes to the 
temperature spectrum, and place limits on the amplitude 
of power coming from cosmic strings. Small-scale tem- 
perature data from ACT constrain the density parameter 
associated with early dark energy models to be less than 
3%. 
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APPENDIX 
A. ACT DATA CONSISTENCY 



In order to c heck consistency of the current best-fit cosmological model with previous ACT constraints presented in 
Dunkley et al. ( |2011J , it is useful to consider the plot derived from the best-fit model. This is shown in the left-hand 



22 



panel of Figure |18| The reduced error bars at larger multipoles yield tighter constraints on the positions of the fifth 
and sixth acoustic peaks, providing tighter constraints on the baryon density, ^Ibh^, relative to the WMAP7-only 
constraints. This comparison is shown directly in Figure |3] 
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Fig. 18. — The best-fit spectra from the ACT data. Left panel: the soHd curve indicates the best-fit Ci from the 3-year ACT data set 
(top) and difference from t he best-fit model (bot tom), while the dashed line shows the curves from the best-fit parameters for the ACT 
1-year data, as presented in |Dunkley et al.| l |20TT} . Right panel: the best-fit power spectra for the WMAP7-I-ACT-E and WMAP7+ACT-S 
data sets. The bands around each best-ht spectrum are all models which lie within la of the best-fit. The spectra are consistent at the 4% 
level. At £ = 2500, 4% corresponds to about 4 /xK^, which is about 25% larger than our error bar at ^ = 2500. 



The marg inalized one-dimension al likelihoods for the ACDM parameters are shown in Figure 19 The ~ 0.7a shift in 



9 a from the 'Dunkley et al. (2011) release relative to the 3-year ACT data constr aints presente ahe re can be explained 
by the similar shift m parameter s of the first 1-year d ata maps used in Dunkley et al. ( 2011[ ) and |Das et al. (2011b I 
relative to the maps released with Diinner et al. (2013). This shift was investigated thoroughly, and found to be caused 
by improved noise treatment and map-making pipeline. This shift in 0a caused a negligible (~ 0.03cr) shift in the 
underlying cosmological parameter fl^. 

The best-fit curves for the ACT-E and ACT-S spectra in combination with WMAP7 data are shown in the right- 
hand panel of Figure 18 while the marginalized likelihoods are also shown in Figure ml The ACT-E data on their 

caTa: 



own (in combination with WMAP7 data) prefer a higher value for ftbh^, 0a and the scalar spectral index. Given the 
~ 4% differences in the best-fit Ci curves, one might ask how well the combined fit compares to the fits obtained when 
considering only the ACT-E and ACT-S subsets individually in combination with WMAP7 data. 

The -2 In if value of the best fit to ACT-S-hWMAP7 data is 444.6 relative to 446.6 when combining ACT-S 
and ACT-E, for 480 data points and 10 fitted parameters, leading to 470 degrees of freedom. When considering 
WMAP7+ACT-E, the -21n^ value is 227.7 compared to 229.0 when using the ACT 3-year data set in combination 
with WMAP7. The relative difference in — 21n^ between the south-only and equatorial-only parameters and the 
respective — 21n^ value for the combined fit is approximately 2; hence one might expect a roughly la shift in 
cosmological parameters when considering the two data sets independently. 
One can go further and test for consistency between individual seasons without including data from WMAP, which 
The roughly la shifts in the value of the baryon density and ratio of the acoustic horizon to 



we show in Figure 20 



the angular diameter distance at decoupling can be understood from the fact that the data are taken from different 
patches on the sky with different noise properties. The two parameter sets are statistically independent, and so ACT's 
internal consistency can be quantitatively checked by approximating the likelihood as Gaussian in the cosmological 
parameters. The covariance of the south and equatorial parameter differences is then simply the sum of the south and 
equatorial covariances, as reported by the MCMC chains. We excise ^ from this test as it is clearly prior-driven. For 
the remaining 17 parameters, we find — 8.78 with a probability to exceed this value (PTE) of 0.88. The kinetic and 
thermal SZ amplitude distributions are meaningfully affected by the constraint that they be non-negative; excising 
these two parameters in addition has only a minimal impact, with — 8.14 for 15 degrees of freedom (PTE=0.84). We 
also find that no particular parameter combinations stand out as being inconsistent between the south and equatorial 
data. 

We divide each parameter difference by its standard deviation and appropriately rescale the combined covariance 
matrix (which now is unity along the diagonal); this operation leaves unchanged. We rotate into the eigenvector 
space of the re-scaled covariance and find that the eigenmode distribution is consistent with a Gaussian distribution; in 
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Fig. 19. — Marginalized 1-dimensional likelihoods for the ACDM cosmological parameters, computed for the ACT-E data in comb ination 
with WMAP7, the ACT-S+WMAP7 combination, and the constraints computed using only the AC T l-year d ata, as presented in Dunk ley] 
|et al.| l |20li| |. The dot-dashed curves show the parameters from the ACT 1-year analysis (Dunklcy et al.|2013^ , while the solid and dashed 
hnes are the parameters from the WMAP7-I- ACT-E and WMAP7-I-ACT- S co nstraints respectively. While there are small shifts between 
the two parameter sets, the best-fit spectra, shown by the curves in Figure [Ts] agree to the 4% level, and the parameter sets are consistent 
when considering the full covariance between parameters. 



particular for both the 15- and 17-parameter case, the most discrepant mode falls within the 68% confidence interval 
expected from purely Gaussian statistics. We note that the posterior parameter distributions need not be Gaussian, 
but that Gaussianity does appear to be a good approximation, and that the cosmological parameters derived separately 
from the AGT south and equatorial regions are fully consistent with each other. 

B. GALACTIC FOREGROUND TREATMENT 



As discussed in Dunkley et al. (20131, we marginalize over residual Galactic dust allowing for two different levels 



in ACT-S and ACl'-li; and imposing Gaussian priors on the dust a mplitudes a^e — 0.8 ± 0.2; — 0.4 ± 0.2, after 

We tested that using a bigger or smaller 
out in the equatorial strip a "seagull" -like 



masking bright dust sources in map space as reported in Das et al. 
mask has a negligible effe ct on all the par ameters 



structure (see Figure 7 in Das et al.||2013 ) 
th e equatoria l map. Finally we cross corre 



Moreover we 

and again found that p arameters do not move if we leave this region in 



Das et al. i2013 



ate ACT and IRIS maps ( Miville-Deschenes fc Lagache| 2005), as described 
and subtract a dust template (rescaled by the cross-correlation coethcient derived through this 
cross-correlarion) to our maps, before computing the power spectrum. We compared parameters obtained with or 
without the dust template subtraction and found the impact on parameters to be negligible. 

C. CALIBRATION AND BEAM UNCERTAINTY 

The calibration factors of the 148 GHz and 218 GHz ACT-S and ACT-E spectra are allowed to vary in the cosmo- 
logical analysis, yielding a 1 (2)% calibration error for the 148 (218) GHz spectra respectively. This error is similar, 
or small er than the calib ration error obtained from calibrating the 148 GHz maps by cross-correlating with WMAP7 
maps (j Hajian et al.|2012 ), which gives an error of 2%. We apply a prior on the calibration parameters for the 148 GHz 
and 218 GHz spectra ot 1.00 ± 0.02. We test whether removing the prior has an effect on the parameters. This shift 
is negligible for the cosmological models considered here. 

Similarly we test for a dependence of the cosmological parameters on the beam error as measured by planet obser- 
vations, by doubling the assumed beam error and comparing this to the case where the beam error is ignored. We 
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Fig. 20. — Marginalized 1-dimensional likelihoods for the ACDM cosmological parameters, computed for the ACT-E data (without 
WMAP7), and the ACT-S data. The spectral index Us and optical depth r are held fixed, and all calibration parameters are set to their 
best-fit values from the ACT-E+WMAP7, ACT-S+WMAP7 runs respectively. We find that the data sets are statistically consistent when 
considering the full covariance between the parameters. 




Fig. 21. — ACDM parameter consistency between the full likelihood and the marginalized Ce bandpower likelihood. The distributions 
are consistent between the two likelihood treatments. 
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find that doubling the beam error shifts the peaks in the foreground parameters by < O.lfi in all cases, relative to the 
constraints when no beam error is applied. 

D. MARGINALIZED CMB LIKELIHOOD 

To compare results to the full mu lti-frequency likelihood , we cross-check constraints on various models using the 
simpler 'CMB-only' likelihood from Dunkley at al. (2013), which first marginalizes over the secondary parameters 
to generate CMB band powers and a covariance matrix. This marginalized likelihood is then used to constrain 
cosmological models directly without requiring any further sampling of nuisance parameters. 




Fig. 22. — One-dimensional likelihoods for the ext ension parameters in various models. The two likelihoods yield consistent results. The 
solid and dashed lines are the same as in Figure [2T] 



D.l. ACDM model comparison 



In Figure 12 of Dunkley et al. (2013) we show that the secondary parameters estimated using this marginalized 
foreground likelihood and assuming the ACDM model are consistent with those obtained from the full multi-frequency 
case. There is a slight (Ict) shift in estimated kSZ amplitude, which shifts due to its correlation with the primary 
CMB, however the constraints are still consistent with an upper limit on the kSZ amplitude. Figure [21] shows the 
marginalized one-dimensional likelihoods for parameters in the ACDM cosmological model. The shift oetween the 
constraints from the full likelihood and those from the 'CMB-only' likelihood are negligible. 

D.2. Extended model comparison 

In addition to the test for the ACDM model, we also test for the consistency between the full likelihood and the 
marginalized CMB-only likelihood for the extended models including the number of effective relativistic degrees of 
freedom, Ncs, allowing for an estimated lensing amplitude, A^, including running of the spectral index, dn g /d In fc 



and variation in the fine- structure constant a/ao- The one-dimensional likelihoods are shown in Figure 22 where 
the maximum likelihood values are consistent between methods; the shifts in the mean values are less then 0.05cr 
in all cases. Both likelihood codes are available on LAMBDA; these tests indicate that the CMB-only likelihood 
returns robust results compared to the full likelihood, and may be used in cases where only the primary cosmological 
parameters are of interest. 



